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ALTERNATIVE ORTHOGONAL RATIONAL FUNCTIONS 
ON A HALF LINE 

VLADIMIR S. CHELYSHKOV * 


Abstract. System of alternatively orthogonalized rational functions of Jacobi type on the half line [1, oo) is 
defined and its properties are established. Three subsystems of proper and mixed systems of rational functions with 
nice properties are presented. 
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1. Introduction. Interest in approximation on a semi-infinite interval by algebraically 
mapped classical orthogonal polynomials on a finite interval originates from 0,fD, and such 
an approach was later developed by many authors (see, for example ( 6 ), 171). 

By following this path, we present system of alternative orthogonal rational functions of 
Jacobi type defined on a half line. We find that the functions possess all the same properties 
as the alternative orthogonal polynomials do. Also, we show that algorithmic capabilities of 
approximation by constructed (both proper and mixed) subsystems of rational functions are 
equivalent to those of the classical orthogonal polynomials. 

2. Bidirectional orthogonalization and the Jacobi polynomials. Procedure of bidi¬ 
rectional orthogonalization of a sequence of functions and alternative orthogonal polynomials 
were introduced and studied in 0 , 0 , 0 . 

Following 0 , let us consider two systems polynomials 


f'r>( a ’P) 

t' nk 


(*)}L 


a > — 2 , 


P>~ 1, 


and { P ifc ,/3) (£)}/-„, > “I 


( 2 . 1 ) 


that are orthogonal with respect to the weight function w(x;a,(3) = x a {l — x)& on the 
interval [0,1]. The systems are constructed by applying inverse and direct orthogonalization 
algorithms beginning with k = n. The polynomials in ED are related by identities 


' nk 


M = p"A(x-') 
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with 0 < k < n and 7 = —a — fi — 2 . 


Also, 


and 


V^\ x )=x k P^ + _T m ^)> 0 <fc< 


Pi a /\x) = x n P^\x) k> 


(2.3) 


(2.4) 


where P 0 m /3) (x) are the Jacobi polynomials on the interval [0,1]. 

Identities \2.2\ - \2A\ were used in 0 to derive classical properties of the alternative 
Jacobi polynomials P^j.' {x). 
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3. Definitions and properties. We obtain alternative orthogonal rational functions on 
the interval [1, oo) making use of key identities (2.2 1 and other properties of the polynomials 


(«. 0 ) 


nk 


X 


) and (a;) established in ||4l . 


{.<*,&)( 

nk * 

First definition and orthogonality. We introduce 

w(x; a, ff) = x a \l — x\ 13 

as a unification formula for the weight functions on the intervals [0,1] and [1, oo) and define 
alternative orthogonal rational functions as follows 

n nk\ X ) = V nk f>] ( X ~ l )> 


0 <k<n, /3 > — 1 and 7 + (3 < 0. 

The functions Tl[^ lj \x) satisfy the orthogonality relation 

/ OO 

(x - 1 fx' y n^ 0) (x)TZ^ ) {x)dx = r^ P) 5 kU 

T(n + k — 7-/3) T(n — k + /3 + 1) 


where 
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’ nk T(n + /c — 7 ) (n — k)\ 

Invariance. For pgN orthogonality ( |3.1| ) holds the invariance 

(7+2p,/3) __ (7.^) 

' n-\-p,k-\-p - nk ’ 


and 


Integral. In addition, one can evaluate 
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Recurrence Relations. With 


(n—k)l T(n — 7 ) 


T^n,n-i{ x ) = “(7 + /? - 2n + 2)a; _ ^ n_1) + (7 - 2n + l)a; _ " 

the functions lZ^^\x) satisfy the three-term recurrence relation 


(3.1) 


(n — k + 1)(7 + /? - n - k + 1)(7 + /? + 2k)TZ ( ^f} 1 (x) 
= (7 + /3 — 2 fc + 1 )[—(7 + (3 - 2k + 2 )(7 + /3 — 2fc)x 
-(7 + /3 - 2 n )(7 — 2k + 1) - 2(n - fe)(n - k + 

-(7 - n - k)(fi + n - ft) (7 + $ - 2k + 2)U^f_l 1 {x). 
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Second definition. Identities (2.2 1 can also be represented as ( cf. 0 ) 

^nk 0 \ x ~ l ) = P-n%(. x ) = xP -^+ll- {k + !)(*). k = n,n—l,..., 0, 
which results in definition 

K-nk^ix) = P-n 0 -k(x). 

We employ ( |3.2| > for deriving results given below. 

Rodrigues’ type formula and integral representation for the rational functions. 


(3.2) 


, t _1 \n-k„n-j(„ _ i \-/3 jn-i 

= ( } d-t)! 


with k = n,n — 1 ,..., 0 , and 


P-nk^ix) = 


(_l)n-fc j.n-7 


g'l—n—kfg _ 2)/3+ n— fc 


27ri (x — \)P J (z — x) n ~ k+1 
c 


dz. 


Here C is a closed contour encircling the point z = x. 

Differential-difference relations. 

(2k + 7 + 2)x(x - P}^Pntf\ x ) = [ n ( n + 7 + P + 2 ) + k(k - /?) 
dx 


—k(2k + 7 + 2)x]7tL^ i .’ / ^(x) + (n + k + 7 + /3 + 2)(ra - k + /3)7£^.+ x (x) 


(7>/3) 


and 

(2k + i)x(x - 1 )~j“^ r nfc^( a: ) = [—n(n + 7 + /3 + 2)-(fc + 7 + /3 - l)(fc + 7 + /? + 1) 


+(2k + 7 )(fc + 7 + l)x]TZ^\x) - (n - k + l)(n + k + 7 - 

Differentiation formula. 
d 


dx 


Pnk 0 ) (x) + -P { nk 0 ) (x) = (n + k + 7 + l)7^f + 1 ) (x). 


Differential equation. The function y(x) = (x) is a solution to the differential 

equation 

x 2 (x — 1 )y" + x((j + /3 + 2)x — 7 — 1 )y' — [k(k — 7 — /? — l)x — n(n — 7 )]y = 0 


with ( 7 , /?) in the domain specified below (see Fig. 6.1 1 . 
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4. Gaussian quadratures and discrete orthogonality. Non-trivial zeros of i( x ) 
provide abscissas of alternative ( 7 , /3)-Gaussian quadratures that are exact for x~i with 


2k — 1 < j < 2n, 0 < k < n. 


(4.1) 


Inequalities (4.1 1 grant discrete orthogonality, but a chosen quadrature is exact for k<j< 
2 n, iff k = 1. We call the quadrature with k = 1 the alternative (7, /3)-Gaussian quadrature 
for rational functions. 

Similarly, one may also find that orthogonal rational functions defined by direct algo¬ 
rithm of orthogonalization as R^ l3 \x) = generate Gauss-type quadrature with 

the zeros of R^k-n'^ (p)- The quadrature is exact for x J with 

2 n < j < 2k — 1, k > n. 

It becomes the (7, /'Ij-Gaussian quadrature for rational functions iff n = 0. 

5. Proper alternative orthogonal rational functions. Let n € N and 

n^ ] (x) = {R^\x)}U n 

be a system of proper rational functions. This orthogonal system is of special interest inas¬ 
much as its terms are quadratically integrable functions: 


/ o° >2 

dx<oo. 


(5.1) 


The associated orthogonal function 7 (x) does not satisfy (5.1 1 , but the zeros of (x) 

are abscissas of the alternative (a, /3)-Gaussian quadrature for proper rational functions. Ab¬ 
scissas and weights of the quadrature can be recalculated from the standard (a, /3)-Gauss 
quadrature over a closed interval. We also find that TV^ 3 ' 1 (x) possesses discrete orthogonal¬ 
ity property with the abscissas chosen. 

All of this makes the set 


n 


(7,/3) 


UK. 


h,P) 

nO 


a system of its own kind. 

6. Singular orthogonality and marginal systems. Let 



7'G {7 | -l<7 + /3<0, —1 </?}. 


(6.1) 







ORTHOGONAL RATIONAL FUNCTIONS 


5 


In this exceptional case the improper rational function lZ^ 0 ’^\x) is not square-integrable 
with respect to the weight function w(x\ 7 '. B ) and therefore non-normalizable. However, it 
is consistent with orthogonality relations (3.1 1 for l > 0. Thus, 1Z^ 0 ’^\x) is a singular term 
of the set 


(x) 


which we call a marginal system. Below we describe two noteworthy orthogonal systems of 
such a kind. 


7. A-kind rational functions. Considering (6.1 1 we find that functions lZ^ 0 (x) = 


K { n0 m (x) are the shifted Legendre polynomials of the reciprocal 1/x, and we count TZ^ 0 (x) 
as the associated term for the system of proper rational functions 

K(x) = {ii£ k (x)}U n , n£ k {x) = n { - k lfi \x). 

The functions TZ„ k (®) obey the orthogonality relations 


^nk 

l^nk(x)n^ k {x)dx = 


k = 0 , 1 , n, 1 = 1,2, ...,n 



and 


f°° 1 1 

J -T^nk(x)dx =-, k = 1 ,..., n. 


They can be calculated using the expansion 

n—k 


K\ 


M = t (-D 5 ("/)(" t-V ) 1 = 0.1. 


or the three-term recurrence relation 

n£ n (x) = x~ n , <„_:(*) = ( 2 n - l)x~( n ~V - 2 nx~ n , 

(2k + 1 )(n + k)(n — k + 

= 2k[(2k - 1 )( 2 k + l)a; - 2 (n 2 + k 2 + n)\R^ k (x) 
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-(2 k - 1 )(n- k)(n + k + l)K^ k+1 (x). 

We also find that y(x ) = 7t// k ( x ) a solution to the differential equation 

x 2 (x — 1 )y" + x 2 y' — [k 2 x — n(n + l)]y = 0. 

The set 

(t^U^q) (x) 

provides regular Gauss-Rational quadrature over a half line, which is exact for for x _J , 1 < 
j < 2 n with the weight 1/x, and discrete orthogonality for 'R.^(x). 

8. T-kind rational functions. Let notation {a) n represents falling factorial. Consider¬ 
ing ( |6.1[ ) we find that 

^noO) = n -/( n ~ h)n n nO~ 1/2) ( X ) 

are the shifted Chebyshev polynomials of the reciprocal 1/x, and we count 7£j 0 (x) as the 
associated term for the system of proper rational functions 

'ft-nO) = {K-Z k { x )}U n 

T^Zk( x ) = c nkR-Zf 1,2 \ x ), c nk = (n - k)\/(n -k - l/2)„_ fc , 

n £ N. The functions IZ^k ( x ) obey the orthogonality relation 



°° K k {x)nl h [x) = 1 (2n + k + 1 — 1)!! (2n-fc-Q!! 

y/x — 1 2 (k + l) — 1 (2 n + k + l — 2)!! (2 n — k — l — 1)!! U 

fc = 0,l,...,n, Z = 1,2,n, 



(2jfe-3)!! 
(2 k)\\ 


2m r, 


k = 1, 2 ...n. 


and 
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They can be calculated using the three-term recurrence relation 

nl n {x) = x ~ n , nl^ix) = (4n - 3)at" (n_1) - (4n - 2)®-", 

[2(?r + fc) — 1] [2(n — fe) + 1] (4fc + (x) 

= (4k - l)[(4fc - 3)(4fc + 1)® - 2(4n 2 + 4fc 2 -2k- 1 )]K% k (x) 

-4 (n - k)(n + k)(4k - 3)TZj lk+1 ( x ). 

We also find that y(x) = H^ k (x) is a solution to the differential equation 

x 2 (x — 1 )y" + x(\x — 1 )y' — [k(k — \)x — n 2 ]y = 0. 

The set 

U 7^-no) ( x ) 

provides Chebyshev-Gauss-Rational quadrature over a half line and discrete orthogonality 
for 'R.l(x). 

In addition, a fast Fourier transform can be employed for computations. 

9. Mixed systems of alternative orthogonal rational functions. Let 

7 + /3<-l, /3 > —1 


and 


iz 


h>P) 


( x ) = {^nk^ ( x )}k=n 


be a mixed system, i.e., the system that contains an improper rational function with k = 0. 
Such a system may be of interest for approximation of functions that have a finite value at 
infinity. 

We introduce alternative Legendre rational functions 

n n (x) = n { ~ m (x), 


considering them as the simplest subsystem of TZ 


(7>/3) 


(x) with nice properties. Functions 


Tln(x) = {7 Z nk (x)}° k=n 


obey the orthogonality relation 

f \lZnk(x)TZ nk (x)dx = Sk ' 1 , k,l = 0,1, 

Ji x z k + l + l 

and their properties directly follow from the properties of alternative Legendre polynomials 
described in 0. The system 1Z n (x) provides Radau-Rational quadrature over a half line with 
one of the abscissas at infinity; also, it possesses discrete orthogonality for {7 Z nk (z)} Un- 
Abscissas and weights of the quadrature can be recalculated from the standard Radau quadra¬ 
ture over a closed interval. 
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10. Conclusion. Orthogonal functions presented in this preprint preserve hierarchy of 
algebraic decay of the sequence xr :l , j = n, n — 1..... 1 at infinity. This property may be 
considered as an advantage for approximation on a half line, in particular if the asymptotic of 
a function is not known exactly. Weak formulation for problems involving differentiation may 
be required. Also, the systems may be of interest for studying energy cascades in nonlinear 
problems. 
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